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Abstract 

We demonstrate an exact local transformation which maps a purely Fermionic manybody system 
to a system of spinfull Bosons and spinless Fermions, demonstrating a possible path to a non-Fermi 

liquid state. We apply this to the half-filled Hubbard model and show how the transformation 

■ 

maps the ordinary spin half Fermionic degrees of freedom exactly and without introducing Hilbert 

O ■ 

space constraints to a charge-like "quasicharge" fermion and a spin-like "quasispin" Boson while 



preserving all the symmetries of the model. We present approximate solutions with localized charge 
which emerge naturally from the Hubbard model in this form. Our results strongly suggest that 
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| charge tends to remain localized for large values of the Hubbard U. 
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The Hubbard model is a paradigm for one of the most difficult quantum manybody 
systems; strong interactions competing with the kinetic energy in a dense electron gas. These 
problems must be explored either analytically by relatively poorly controlled approximate 
methods or by methods such as Monte-Carlo, dynamical mean field theory (DMFT) or finite 
size calculations. A major difficulty is the absence of a framework for such systems that 
competes in simplicity and rigor with conventional Landau Fermi-liquid theory. 

A particularly intriguing scenario for non-Fermi liquid behavior is "spin-charge sepa- 
ration" with a low-energy effective theory where spin and charge excitations are decou- 
pled. This idea has attracted widespread attention because it migh provide a route to 
high-temperature superconductivity in a system with strong local repulsion. In 1-D it is 
known through the use of "bosonization" that spin-charge separatation does occur at low 
energy.^] In higher dimensions, much less is known. The usual approach to spin charge sep- 
aration is to construct electron operators as composite objects consisting of a charge part, 
the "holon", and a spin part, the "spinon". To avoid unphysical states contraints have to 
be included which are then treated by approximate methods. Q] 

Here we shall construct electron operators as composites of charge-like and spin-like 
operators that do not give rise to any unphysical states and thus avoids introducing any 
constraints. An exact local mapping which is valid in all dimensions will be shown to exist 
from a purely Fermionic manybody system to an equivalent system of Fermions and spin- 
like operators obeying Bosonic commutation relations between different lattice sites. We 
shall explore this with the 2-D Hubbard model, for which the resulting mapping bears a 
resemblance to the mappings of the local Hilbert space achieved by introducing spinons and 
holons. 



Transformation to quasiparticle operators. 

Let n r be the particle number at site r, and n-f and be respectively the number of 
up and down electrons: n r = n^ r + n^ r . Since the Hubbard interaction obeys n^ r n^ r = 
|(n r — 1 + (n r — l) 2 ), by choosing the chemical potential appropriately we can replace the 
usual Hubbard interaction with the particle-hole symmetric interaction ^ U {n r — l) 2 . We 
thus write the Hubbard model in the particle-hole symmetric form H = —t J2{rr')( c t,r c s r ' + 
CC) + U ^ r | (n r — l) 2 where the sum is over all sites and nearest neighbor bonds counted 
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once. 

We now demonstrate an exact transformation of the local operator algebra of the Hubbard 
model to new Fermionic CP 1 "quasicharge" c r and SU (2) "quasispin" operators q\ that obey, 
respectively, Fermi and Bose statistics. The operators are given by 

K = c{,r(! - n i,r) + (-l) r C Tjr n ijr (1) 

at = (c[ r -(-iyc^ r )c Ur (2) 

Qr = (?r + ) t 
Qr 2 n lt r 

where — V is ±1 depending on which sublattice site r is on. We define the x and y com- 
ponent of quasispin by qf = | (g+ + q~) and q% = ^ (g+ — q~). The factor appearing as 
(— l) r can in fact be chosen as an arbitrary phase, but is fixed here to preserve symmetries 
specific to the Hubbard model. This "nonlinear" transformation generalizes earlier work on 
transformations that were strictly canonical, and hence did not alter commutation relations 
of the fermi operators. jjj 0] 

The new operators obey the following algebra, which can be verified by straightforward 
manipulations: {c r ,c\,} = 5 r y , {4,4'} = > [4>#'] = > [q^qI'] = i5 rr> Y,k e ijkQr 

We thus find that q l r1 c\ and c r are respectively independent spin half bosonic and "spin- 
less" Fermionic operators. The usual operators are given by 



c[ r = c r {\ + q z r) + (-l) r 4(3 - 9?) ( 3 ) 
c{ r = q~ (c r - (-l) r cj). 

We define the quasicharge operator as n{. = c\c r . We find the following relations with the 
usual operators expressed in terms of ordinary Fermions: n r = 1 — 2n{.q*, s\ = (1 — n{. ) q\ 
and n{ = (n r — l) 2 where s\ = f52 a 8 c oi,r a a/3 c p,ri w ith a % the Pauli matrices. We also 
define the local "pseudospin" operators p l r = n{q\ which are the generators of the SU(2) 
algebra which correponds to "rotations" between the empty and doubly occupied states. |5( 
The formal resemblance to spin symmetry has led to the name pseudospin. The total z- 
component of pseudospin can therefore be seen to be half the number of doubly occupied 
sites minus the number of empty sites, which is precisely the charge relative to half filling. 
Pseudospin is a symmetry of the particle-hole symmetric Hubbard model but is broken by 
a chemical potential. 



We note that, aside from phase factors, the quasicharge operators always create even 
charged states from odd and vice versa, as they must since they are Fermionic. However, 
there is a precise phase relationship which is not so easily interpreted that is enforced by 
the commutation relations and symmetries. 

Our Hubbard model can now be rewritten exactly in terms of the quasiparticle operators 

as 

H = t(T + T 1 +T. 1 ) + Uh u (4) 

with hu = | J2 r $rK an d 

To = \ (1 + 4<Zr • q r >)(ctc r , + CC) (5) 
Ti = \ E "l r (1-4^-^0(44,) 

< r,r' ) 

and T_i = T\. In this representation, the symmetry under the entire SU{2) quasispin rota- 
tion is manifest and we note that the Hubbard interaction becomes simply the quasicharge 
number operator. 

We define total quasispin Q , spin S and pseudospin P respectively by Q = J2 r Qr 
etc. P, Q and S all obey the SU(2) algebra [Pi,Pj] = iJ2k e ijkPk etc. We see from the 
definitions above Eq. HJthe relationship, p l r = n{.q l r and s l r = (1 — n{. ) q l r . Thus Q = P + S, 
i.e. quasispin can be exactly split into pseudospin and spin, with nf. and (1 — n() providing 
the projection operator of Q into either spin or pseudospin. Since Q commutes with c\ 
and the Hubbard Hamiltonian H depends on q r only through rotationally invariant terms, 
we conclude that [H,Q] =0. We further know that the Hubbard model commutes with 
ordinary spin whereby [H, S ] = 0. We can therefore conclude that [H, P ] = 0, which 
confirms the well known invariance under pseudospin rotations. 

The operators T±io are identical to the operators defined in Ref. 

7] 

and used in a 

perturbation expansion in t/U. It is clear from their definitions that T±i couples different 
Hubbard bands defined by the sectors given by different values of (hu) whereas To does 
not. We observe [T 0i ±i , Q] = 0. It can be inferred that [T ±i , P] = since T ±i is known 
to commute with ordinary spin. 

We now define the usual canonical transformation S which splits the Hamiltonian into 
pieces that preserve total quasicharge. We define S = i(T_ 1 — Ti) + i [ T , (Ti + T_x) ] 
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and it straightforward to show that [S, hu] = (T_i + Ti), whereby it follows that H e ff = 
e us / u H hub e- us / u = H hub +it[S,H hub }/U- t 2 [S, [S, H hub }]/(2U 2 ) + ... = H q£ + H q + 0(t 3 /U 2 ) 
splits into two terms. The first term H q is the ordinary spin interaction in the "spin-only" 
sector 

H q = jJ2(Qrqr'~\) (6) 

(rr'} 

where J = 4t 2 /U and the second term H qc couples quasispin and quasicharge 

t 2 

Hqc = Yl ( UflU + tT + ^j( nf r+ n f r>)(qrq r > ~ \)) + 

( rr' } 

t 2 

+ 177 ^A" + d t" d r) (! - <?r<?r' + q r q r " + 



4U 

( rr'r" ) 



- q r q r " + zq r • (q r i x qr r ») .) (7) 



The last summation runs over all three-site neighbors connected by links. Terms of order 
tfjU which are already present to order t or U are ignored. According to the results of Ref. 
pi, to all orders in t/U the operator S which diagonalizes total quasicharge is a polynomial 
of T 0j -i-i, and hence to all orders in perturbation theory, the transformation by e lts l v will 
preserve all the symmetries Q and P in this formulation. 



The spin-only sector 

Under this canonical transformation the eigenstates of H e ff will map to the eigenstates 
of the Hubbard model represented by the bare quasiparticle operators through the transfor- 
mation e lts l u '. We let | ^ e ff) be eigenstates of H e ff and | ^ ) = e tts / u \ ^ e //)- Adapting 
Ref. [3| to our problem, there will be a "zeroth Hubbard band" or "spin-only sector" of 
solutions corresponding to n{ \ ^ e ff ) being identically zero so that H e ff reduces exactly to 
H q . We conclude (H q £)^ ^ is identically zero and the effective Hamiltonian of the zeroth 
subband to order t 2 /U is exactly the Heisenberg model. 

What can we conclude from symmetries about the properties of the zeroth Hubbard band? 
It must be remembered that N[ ot = ^~J r n? is not a "good quantum number" i.e. is neither 
invariant under the symmetries nor under e tS . We have therefore no information about 
(n{,)^ from symmetry arguments alone. However, we see that nl \ ^ e //) — throughout 
the entire lattice in the entire zeroth Hubbard subband. We can therefore conclude that 
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(P)y = 0. Since P commutes with S, this assertion survives to all orders in perturbation 
theory and we can conclude that (P)q, = 0. In fact quasispin and ordinary spin coincide. 
Since ( n ) = ( 1 — 2P Z ) is equal to physical charge, we see that the spin-only sector has 
indeed charge density corresponding to half filling. 

Extended states in the charge-only sector 

The pseudospin is given by P 2 = ^2 rr , q r • Qr 1 n/^v If there is exactly one quasicharge, 
we conclude from the fact that q 2 . = 3/4 that ( P 2 ) = |. We have thus identified the lowest 
nonzero quasicharge sector with the spin half representation of pseudospin, exactly what we 
need to associate it with one added positive or negative physical charge. 

We can derive a set of extended charge states in the dilute limit by a mean field treatment 
of Eq. [3 replacing q r ■ q r r — ( q r ■ q r > ) H . In two dimensions, these have the spectrum 

E k = ~U — 4to(cos k x + cos k y ) + 

t 2 

— ((cos2fc x . + cos2k y )(l - e 2 ) + 

2 cos k x cos k y (1 — e^)) (8) 

where a = |(| + e i)| an d where t\ = (q r q r '), an d e 2 denote nearest neighbor and 
longer range Heisenberg correlations. Using the numerical value e\ = —0.33972(2) [s| we 
find a = .08972. We find a quasicharge spectrum bounded by U/2 ± 8at + 0(t 2 /U) and 
shown in the shaded region of Fig. |21 

The Nagaoka instability and localized states 

The present calculation makes the "Nagaoka theorem" natural. 6] Consider Eq. 0] in 
the case where the quasispins are perfectly aligned along the quasispin z axis. In this case 
q r ■ q r i = i, the operator T±i will be identically zero and the kinetic energy T becomes 
precisely ^2 rr , clc r ,+cl,c r . These are quasicharges moving independently of the ferromagnetic 
quasispins and are ordinary holes in a perfect 2-D ferromagnet which the kinetic energy 
— 2t(cos k x a + cos k y a) + \U in the thermodynamic limit. 

The Nagaoka theorem states that as U — > 00 a single charge causes the ground state 
of H to become ferromagnetic. This Nagaoka ferromagnet with charge one has Hubbard 
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FIG. 1: (a,b) The two smallest quasiferromagnetic clusters is shown. The shaded regions include 
the sites on which charge is approximately confined. 

energy Ef erro = — At, whereas the antiferromagnetic solution with a single quasicharge 
has energy E anti f erro = 2VL.J (ei — |) + U/2 — Sat where Q is the number of lattice sites in 
the system. We therefore find that Ef erro < E anti f erro when the inequality U > Qt(l + 
2a)/(l — 2a) holds. We see that whether or not there is Nagaoka ferromagnetism depends 
on which order we take U and Q to infinity. 

The argument can be used to estimate when the antiferromagnet becomes unstable to 

n 

ferromagnetic ordering. (See also Ref. [9( ). Let us consider a quasicharge density of 
p <C 1. We see that if p obeys p < t/U (1 + 2a) J{1 — 2a) the dilute ferromagnet will be 
energetically favored over a collection of extended quasicharge states, essentially reproducing 
the argument by Ioffe and Larkin.j^] 

Using the previous ideas, we shall study localized states in a 2-D antiferro magnet. As we 
have seen, a quasicharge hops freely on a ferromagnetic bond where q r • q r < = |. In a mean 
field treatment of the quasispins, to linear order in t a quasicharge cannot hop an a bond 
with antialigned quasispins since ( q r ■ q T i ) = — |. 

Starting with a patch of locally Neel ordered quasispins oriented along ±q z we consider 
a quasispin state with nearby spins flipped relative to the Neel state on diagonally adjacent 
sites of the spin-down sublattice. We note that in the presence of a quasispin ferromagnet 
oriented along +z, quasicharge corresponds to physical charge measured from half filling. 

The core of the region will thus be a ferromagnet surrounded by a closed boundary of 
antialigned spins. Outside this patch, quasispins are allowed to relax, so that far away 
the system will be a Heisenberg antiferromagnet. The two smallest such configurations are 
shown in Fig. d 
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We identify the quasispin of the state by constructing the Neel state with the same 
rotational symmetry around the cluster. The smallest cluster (Fig^a) has one flipped spin 
thus has Q = \ whereas the cluster with two flipped spins (Fig. [TJb) has Q = 2. We now 
add rif quasicharges to the cluster. According to the definition of pseudospin, in the case 
where quasicharge is confined to a quasiferromagnetic region with q r ■ q r i = ■? + ~5 rr / , it 
can be shown that P 2 — -y(-y + 1)- A cluster with rif units of quasicharge will thus have 
pseudospin P = ^ n f an d, from the identity Q = S + P, physical spin S = \Q — \rif\. 
The clusters with no quasicharge have pseudospin zero and physical spin equal to quasispin, 
whereas the charge one cluster will have pseudospin P — | and spin S = Q — |. 

The energies of these states are approximate but the spin and pseudospin quantum num- 
bers are exact. As quasicharge is added in the cluster, quasispin, spin and pseudospin con- 
tinue to be good quantum numbers. Since quasispin is conserved upon adding quasicharge 
each such configuration has a well defined value of quasispin independent of physical spin 
or charge and can be calculated in the effective theory. Quasicharge, however, is not a good 
quantum number, and will change upon transforming back the the physical state with e tts l v '. 

We estimate the energy by keeping only the part of the kinetic energy linear in t and ignore 
the spin exchange energy beyond the antiferromagnetic boundary of the cluster. For the 
smallest Q = ~ cluster we see sixteen bonds that are strongly affected by the localized charge. 
Four bonds are ferromagnetic with zero energy and twelve bonds are antiferromagnetic with 
energy J x (— |), i.e. they retain about 85% of the energy of the unaffected Neel state 
J x (ei - J) « -.59J. 

The energy of these states is compactly written with the spectroscopic-like notation ®Eg 
where Q is the quasispin, c = 2P Z is the charge relative to half filling and S is the spin. 
Thus ^-E? is the energy of the the smallest uncharged cluster measured with respect to 
the Heisenberg ground state. We find = J(12 x (-1) + 16 x \e x + \\) w 13.74 t 2 /U. 

The five kinetic energy eigenvalues of this cluster are t(-2, 0, 0, 0, 2). We thus find ^ = 
U/2 + ? El - 2t » U/2 + 13.74 t 2 /U - 2t. For the Q = 2 cluster we find a kinetic energy 
for a single charge — V~Gt; there are 16 antiferromagnetic bonds and 8 ferromagnetic bonds 
giving 2 E\ = U/2 + 2 E\ - V6t » U/2 + 24.61 t 2 /U - 2.451 The energy of these states 

2 

are plotted in Fig. El and compared to the energy of the band of extended states from Eq. 
IH1 We find that for large U the clusters with localized charge are energetically favored over 
the extended states. As U becomes larger, larger clusters lower the energy in accordance 
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FIG. 2: Energy of the extended and the two simplest localized single-charge states, 
with the Nagaoka theorem. 

Conclusions 

We have demonstrated an exact local mapping from a Fermionic many-body system 
to system of interacting Bosons and Fermions that does not introduce constraints. The 
mapping provides a paradigm for a non-Fermi liquid behavior for correlated Fermions valid 
in any dimension. Applied to the 2-D Hubbard model the transformation makes precise 
earlier attempts that relied on approximations and constraints to discuss spin and charge 
as composite operators. Our calculations supports the idea of low energy localized states 
where spin is bound to charge in the limit of large U. 
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